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ABSTRACT 


In  this  report  a  mode-controlled  giant  pulse  laser 
oscillator  and  amplifier  are  described.  The  oscillator 
employs  a  temperature-controlled  resonant  reflector,  a 
dye  Q-switch  laser  and  a  2  mm  aperture  to  assure  single 
mode  operation.  The  result  of  preliminary  threshold 
measurements  on  a  few  different  ruby  samples  are  pre¬ 
sented.  Equations  which  describe  the  propagation  of  the 
gaussian  mode  are  examined  in  order  to  answer  some  impor¬ 
tant  questions  concerning  focusing  inside  dielectric 
media.  The  potential  usefulness  of  the  scanning  elec¬ 
tron  microscope  as  a  tool  for  examining  surface  damage 
is  pointed  out  by  the  presentation  of  a  few  interesting 
photomicrographs.  The  theoretical  effort  is  concen¬ 
trated  on  the  effects  of  strong  optical  fields  on  rela¬ 
tively  low  densities  of  electrons  in  the  conduction  band 
of  a  transparent  polar  crystal.  The  present  studies  have 
led  to  agreement  with  the  general  opinion  that  these 
electrons  are  a  key  link  in  the  optical  pulse  damage 
process.  However,  the  present  work  indicates  that  signi¬ 
ficant  electron  heating  above  ambient  lattice  tempera¬ 
tures  does  not  occur  at  the  optical  intensities  in  ques¬ 
tion.  An  entirely  new  role  of  the  electrons  in  laser 
induced  damage  is  suggested. 
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I. 


EXPERIMENTAL  STUDIES  ON  OPTICAL  DAMAGE 


A.  Introduction  and  Summary  of  Results 


The  experimental  effort  during  the  reporting  period  has 
been  devoted  mainly  to  working  toward  obtaining  reliable  single 
mode  operation  of  our  ruby  oscillator  and  sufficiently  high 
power  to  exceed  damage  threshold  in  ruby  and  sapphire  under 
mild  focusing  conditions.  To  these  ends  we  have  carried  out 
a  number  of  refinements  on  a  giant  pulse  ruby  laser  which  is 
essentially  a  copy  of  one  described  previously,1  and  have 
built  an  adequate  ruby  laser  amplifier.  The  main  modification 
in  the  oscillator  is  the  installation  of  a  temperature  con¬ 
trolled  resonant  reflector  which  will  give  reliable  longitudi¬ 
nal  mode  control.  The  resonant  reflector  used  previously  did 
not  employ  temperature  control,  and  thermal  drifting  resulted 
in  poor  reliability  in  obtaining  a  single  longitudinal  mode. 
Calibration  of  the  power  monitoring  detector  has  been  accom¬ 
plished  by  two  independent  means.  Various  time  delay  circuitry 
for  the  oscillator-amplifier  combination  has  been  built,  in¬ 
stalled,  and  demonstrated  to  perform  reliably.  Preliminary 
threshold  measurements  have  been  made  on  a  few  different  samples 
so  that  some  important  comparisons  could  be  made.  We  have  ex¬ 
amined  equations  which  describe  the  propagation  of  the  gaussian 
mode  in  order  to  answer  some  important  questions  concerning 
focusing  inside  dielectric  media.  The  potential  usefulness 
of  the  scanning  electron  microscope  as  a  tool  for  examining 
surface  damage  is  pointed  out  by  the  presentation  of  a  few 
interesting  photomicrographs.  With  the  exception  of  the  tem¬ 
perature  controlled  resonant  reflector,  which  it  is  hoped 
will  be  available  soon,  we  are  in  an  excellent  position  to 
begin  a  variety  of  meaningful  experiments. 


B.  Experimental  Apparatus 


The  experimental  setup  is  shown  in  Figs.  1  and  2.  The 
laser  head  is  shown  in  Fig.  3.  The  oscillator  employs  a  4  in. 
long  x  1/4  in.  diameter  ruby  which  is  pumped  by  two  linear 
lamps  in  a  double  elliptical  pump  cavity.  The  ruby  crystal 
is  water-cooled  by  a  closed-cycle  refrigeration  system  main¬ 
tained  at  0°C.  The  high  reflectivity  mirror  is  coated  with 
a  99+%  reflectivity  high  field  damage  coating  from  Perkin 
Elmer  Corporation.  Q-switching  is  accomplished  with  a  rotating 


1 


2 


cube  prism  placed  such  that  the  beam  enters  and  leaves  at 
Brewster's  angle  when  the  system  is  in  aligrment.  The  tenqper- 
ature  controlled  resonant  reflector  is  desc  ibed  in  more  de¬ 
tail  in  the  following  section.  Collimated  light  from  the 
He-Ne  alignment  laser  is  used  in  conjunction  with  the  auto- 
collimating  telescope  to  align  the  various  elements  in  the 
laser. 

Portions  of  the  laser  beam  are  split  off  in  various  ways 
(see  Pig.  1),  so  that  the  power  output,  near-  and  far-field 
patterns,  and  Pabry-Perot  patterns  can  be  monitored  for  each 
shot.  This  is  accomplished  in  the  following  way.  Light  re¬ 
flecting  from  wedged  beamsplitter  gives  two  diverging 

beams;  one  of  them  hits  the  magnesium  oxide  diffuser,  where 
the  scattered  light  is  monitored  by  a  biplanar  photodiode 
used  as  our  power  monitor.  The  second  beam  from  hits 

ground  glass  screen  G,  where  it  is  photographed  through  lens 
L  and  the  1  m  focal  length  camera  focused  at  infinity. 

This  gives  a  magnified  (^  Sx)  near-field  picture.  Another 
portion  of  the  light  is  removed  by  beamsplitter  Bi  and  hits 
mirror  t^,  which  can  be  placed  in  or  out  of  position  depend¬ 
ing  on  the  use  of  the  alignment  laser.  From  II2  the  light 
either  goes  to  the  Fabry-Perot  interferometer  or  can  be 
partially  reflected  from  wedged  beamsplitter  W2,  where  it 
results  in  a  pair  of  far-field  patterns.  An  example  of  the 
data  is  shown  in  Fig.  4.  A  0.6  neutral  density  filter  is 
placed  near  the  focal  plane  of  the  camera  so  that  the  far- 
field  pattern  and  the  Fabry-Perot  pattern  can  be  seen  at  two 
different  exposures.  The  two  Gian  prisms  are  used  as  a  vari¬ 
able  attenuator  after  the  amplifier.  Beamsplitter  B,  sam¬ 
ples  the  light  to  the  photodiode  which  monitors  the  power 
incident  upon  the  focusing  lens,  which  was  designed  for  mini¬ 
mum  spherical  aberration  (Special  Optics) .  A  photograph  of 
the  amplifier  is  shown  in  Fig.  5. 

The  water  cooled  amplifier  ruby  is  6  in.  long  x  0.5  in. 
diameter,  with  one  end  wedged  relative  to  the  other  by  about 
0.5°.  The  input  end  of  the  amplifier  rod  is  antireflection 
coated  to  minimize  the  chances  of  oscillation  within  the 
amplifier  itself.  The  ruby  rod  is  closely  coupled  to  a  heli¬ 
cal  flashlamp,  which  is  pumped  with  a  power  supply  capable 
of  delivering  8  kJ  in  a  3  msec  pulse.  The  power  simply 
employs  a  pulse  shaping  network  of  20  sections,  each  section 
pumping  for  150  usee. 
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Pig.  4. 

Typical  photo  monitoring  Fabry-Perot 
inter farogram,  near- field,  and  far- 
field  beam  patterns.  The  difference 
in  optical  density  between  the  two 
halves  of  the  picture  is  N.D.  0.6. 
The  free  spectral  range  of  the  inter 
ferometer  is  1.6  cm-1 . 


C .  Longitudinal  Mode  Selection 


A  multiple  element  resonant  reflector  has  been  designed 
to  enhance  longitudinal  mode  selection.  The  optical  resonator 
and  associated  equipment  for  accurate  temperature  regulation 
of  the  system  is  shown  in  Fig.  6.  The  resonator  consists  of 
two  quartz  plates  optically  contacted  to  a  quartz  spacer  with 
dimensions  of  5.14,  3.42  and  3.00  mm,  respectively.  It  is 
well  known  that  except  for  extra  modes  in  the  regions  of  high 
loss,  the  general  effect  of  the  extra  surfaces  added  to  the 
laser  cavity  is  to  impose  a  modulation  on  the  original  modes. 
The  frequency  spacings  of  the  plates  are  2/3  and  1  cm“l. 
Thus,  the  combined  effect  of  the  two  plates  is  to  enhance 
cavity  modes  separated  by  2  cm"^  and  to  discriminate  against 
intermediate  modes,  while  the  surfaces  separated  by  the  spacer 
enhance  modes  separated  by  1-2/3  cm“l.  A  temperature  change 
alters  the  optical  length  of  the  plates,  but  has  little  effect 
on  the  length  controlled  by  the  spacer  because  the  temperature 
coefficient  of  expansion  for  quartz  is  an  order  of  magnitude 
lower  than  the  refractive  index  coefficient.  Thus  the  device 
can  be  thermally  tuned  to  bring  all  of  the  surfaces  into 
resonance  to  provide  a  maximum  reflectivity  of  40%  at  the 
resonant  frequency. 

The  optical  transmission  of  light  from  a  He-Ne  laser 
through  the  resonant  reflector  at  different  temperatures  is 
shown  in  Fig.  7.  The  bandwidth  of  this  laser  is  small  com¬ 
pared  with  that  of  any  of  the  resonator  spacings,  and  its  out¬ 
put  frequency  does  not  necessarily  correspond  to  a  resonant 
frequency  of  the  device;  therefore,  maximum  reflectivity  is 
not  expected.  Measurements  of  this  type  were  made  in  order 
to  gain  information  which  would  simplify  application  of  the 
device  to  the  pulsed  ruby  laser.  The  bandwidth  of  the  ruby 
gain  profile  is  large  enough  to  include  two  resonances  of  the 
reflector  system.  Independent  thermal  control  of  the  ruby 
and  resonant  reflector  will  be  used  to  bring  the  system  into 
resonance  to  provide  optimum  longitudinal  mode  selectivity. 

The  resonant  reflector  is  not  yet  permanently  incorpo¬ 
rated  in  our  system.  This  fact  represents  the  major  obstacle 
to  our  being  completely  ready  to  perform  a  number  of  desired 
damage  threshold  experiments.  There  have  been  a  number  of 
difficulties  associated  with  the  fabrication  and  installation 
of  the  resonant  reflector,  all  of  which  can  be  summed  up  in 
the  result  that  it  will  not  stay  together;  i.e.,  the  optical 
contacting  has  not  been  accomplished  for  a  variety  of  reasons. 
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At  this  writing,  a  set  of  elements  for  the  resonant  reflector 
is  being  reworked  by  the  fabricators  (Crystal  Optics  Research) , 
and  it  is  hoped  that  they  will  be  received,  permanently  opti¬ 
cally  contacted,  within  a  few  days.  If  this  reworking  is 
successful,  we  will  soon  be  able  to  operate  the  laser  as  a 
reliable  single  mode  device. 


D.  Results  of  Threshold  Experiments 


Early  in  the  program  we  had  the  use  of  a  mode  controlled 
giant  pulse  laser  to  carry  out  a  few  preliminary  damage  thresh¬ 
old  experiments.  The  main  purpose  of  these  early  threshold 
experiments  was  to  determine  whether  an  amplifier  would  be 
needed  to  reach  damage  threshold.  We  found,  in  fact,  that 
with  the  oscillator  alone  we  were  not  able  to  damage  the  ruby 
sample  internally  when  we  focused  with  a  50  cm  focal  length 
lens;  however,  a  small  pit  was  formed  at  the  exit  surface  of 
the  sample.  When  a  12-in.  ruby  laser  amplifier  was  used  in 
conjunction  with  the  oscillator,  we  were  able  to  generate  ex¬ 
tensive  bulk  damage.  The  amplifier  was  pumped  to  give  a  gain 
of  about  7  dB.  An  initial  damage  threshold  value  for  a  ruby 
sample  from  Union  Carbide  was  reported  to  be  about  3  GW/ cm2 . 

Having  established  the  need  for  an  amplifier,  we  designed 
and  built  one.  Once  the  amplifier  was  operating  reasonably 
well,  we  were  able  to  obtain  a  rough  comparison  of  damage 
threshold  between  a  typical  Union  Carbide  sample  and  a  ruby 
boule  section  from  Airtron,  which  was  particularly  impressive 
in  its  excellent  optical  and  scattering  properties.  The  Air¬ 
tron  sample  had  been  on  loan  to  C.M.  Stickley  of  Air  Force 
Cambridge  Research  Laboratories.  Airtron  subsequently  con¬ 
sented  to  sell  it  to  Hughes  for  the  damage  study.  The  Airtron 
sample  has  a  damage  threshold  of  about  a  factor  of  two  higher 
than  the  Union  Carbide  sample.  The  measured  values  are 


Airtron  Sample 

Union  Carbide 
Sample 

Threshold  Power  Density 

18.8  GW/cm2 

9.05  GW/cm2 

Threshold  Energy  Density 

609  J/cm2 

295  J/cm2 
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These  measurements  were  made  with  no  attempt  to  maintain 
longitudinal  mode  control,  but  with  good  transverse  mode  con¬ 
trol  as  shown  by  near-  and  far-field  photographs.  (The  far- 
field  beam  divergence  of  the  laser  is  typically  0.7  mrad  (full 
angle),  as  measured  in  the  far-field  photographs;  the  calcu¬ 
lated  diffraction  limited  divergence  for  a  gaussian  beam  of 
1  mm  radius  is  ^  0.5  mrad.  0  -  2X/ttwq.)  The  above  power 

and  energy  densities  were  calculated  using  the  assumption  that 
the  beam  radius  is  the  same  at  the  output  of  the  amplifier  as 
at  the  input  (i.e.,  1  mm).  The  beam  radius  at  the  minimum  was 
taken  to  be  1.1  x  10"2  Cm,  which  is  the  number  calculated  ac¬ 
cording  to  the  method  described  in  the  next  section.  The 
energy  density  was  calculated  using  the  peak  power  times  the 
pulse  width  (32.5  nsec).  The  power  monitors  were  calibrated 
against  a  calibrated  Korad  biplanar  photodiode  and  calorimeter. 


E.  Definition  of  Threshold 


The  above  comparison  of  damage  thresholds  for  the  differ 
ent  samples  must  be  taken  as  tentative  because  the  definition 
of  damage  threshold  is  somewhat  arbitrary.  The  data  which 
led  to  that  particular  set  of  numbers  yield  the  following 
observations: 


Sample* 

Power  Density 
at  Focus, 
GW/cm2 

Remarks 

(4)  Airtron  boule 

7.21 

No  damage 

(2)  Airtron  boule 

7.45 

Surface  pitting 

(5)  Airtron  boule 

18.8 

2  internal  pits 

(3)  Airtron  boule 

36.9 

Extensive  internal 
damage 

(6)  Union  Carbide  Rod 

9.05 

1  internal  pit 

(1)  Union  Carbide  Rod 

39.3 

Extensive  internal 
damage 

The  numbers  in  this  column  refer  to  the  order  of  laser  shots. 
The  sample  was  moved  so  that  a  different  region  was  exposed 
for  each  shot.  After  each  shot  the  sample  was  examined  by 
illuminating  with  a  He-Ne  laser  alogned  to  the  same  direction 
of  propagation  as  the  ruby  laser. 
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Shots  5  and  6  were  taken  to  be  threshold  in  that  a 
very  small  amount  of  internal  damage  is  observed,  but  the 
reproducibility  of  this  observation  is  not  yet  known.  In 
addition,  the  difference  in  power  required  to  give  one  or 
two  internal  pits,  although  believed  to  be  insignificant,  is 
not  really  known.  We  present  these  data  to  point  out  the 
tentative  nature  of  the  reported  threshold  numbers.  It  would 
be  more  realistic  to  state  that  the  damage  threshold  for  the 
Airtron  boule  lies  between  7  and  19  GW/cm 2  and  that  for 
the  Union  Carbide  rod  is  9  GW/cm*  or  less.  We  also  wish 
to  point  out  that  the  surface  damage  at  the  exit  surface  occurs 
at  appreciably  lower  power  than  the  internal  damage;  no  special 
surface  cleaning  techniques  were  used,  however. 

In  a  series  of  experiments  during  the  exploration  of  the 
damaging  capability  of  our  amplifier-oscillator  combination, 
we  placed  a  piece  of  exposed  film  a  few  centimeters  beyond 
the  ruby  sample  to  observe  burn  patterns.  When  extensive  bulk 
damage  occurred,  the  burn  spot  was  about  twice  as  large  as 
that  resulting  when  there  was  no  damage  or  just  surface  dam¬ 
age.  This  probably  results  from  the  increased  scattering  of 
the  laser  light  by  the  damage  sites.  The  effect  of  increasing 
the  angle  of  forward  scattered  light  has  been  reported  previ¬ 
ously.2  The  fact  that  the  damage  is  formed  during  the  giant 
pulse  is  also  substantiated  by  other  workers,1  as  well  as  by 
some  of  our  early  results"  where  we  used  an  STL  camera  to 
photograph  the  damage  in  glass  while  it  was  being  formed. 


F.  Electron  Microscope  Observations 


We  have  used  a  scanning  electron  microscope  to  examine 
the  damage  pits  formed  at  the  exit  surface  of  a  ruby  sample. 
The  particular  sample  is  a  1/4  in.  diameter  x  3  in.  long  rod 
from  Union  Carbide.  Figure  8  shows  several  views  of  the  sur¬ 
face  damage  under  different  magnifications.  The  first  picture 
(Fig.  8(a))  shows  the  rod  under  23x  magnification.  The  dam¬ 
aged  areas  were  caused  by  three  different  shots  of  the  laser.; 
the  one  small  pit  at  the  lower  left  was  caused  by  one  shot, 
the  cluster  at  the  lower  right  by  another,  and  the  large 
craters  at  the  top  by  a  third.  Further  magnification  of  the 
last  two  regions  is  seen  in  the  other  electron  micrographs  in 
the  figure.  Figures  8(c)  and  (d)  show  the  large  crater  at 
the  top  magnified  610x  and  7000x,  respectively.  Figures  8(e) 
and  (f)  show  one  of  the  multiple  pits  on  the  lower  right 
magnified  660x  and  6600x,  respectively.  Electron  micro¬ 
graphs  of  the  undamaged  surface  at  7000x  magnification  show 
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(a)  End  of  ruby  rod 
at  23x. 


Fig.  8.  Scanning  electron  micrographs  of  the  exit  surface 
deunage  at  different  magnifications. 
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Fig.  8  continued. 
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essentially  a  smooth  surface  (Fig.  8(b)).  The  magnified  view 
of  the  large  crater  in  Figs.  8(c)  and  (d)  shows  what  appears 
to  be  a  molten  region  around  the  hole,  whereas  the  views  of 
the  pitted  area  (Figs.  8(e)  and  (f))  appear  to  be  more  like 
multiple  fractures  than  a  melting  phenomenon. 

We  do  not  wish  to  offer  any  interpretation  of  the  photo¬ 
micrographs  at  this  time,  but  merely  present  them  and  point 
out  their  potential  use  as  a  tool  for  studying  laser  induced 
surface  damage. 


G.  Mode  Propagation  Computations 


In  order  to  have  a  more  detailed  understanding  of  the 
behavior  of  the  light  after  it  propagates  away  from  the  laser, 
we  have  calculated  the  beam  size  as  a  function  of  distance 
from  the  laser  and  thence  through  the  lens  to  the  sample. 

Our  interest  in  pursuing  these  considerations  was  to 
answer  the  following  questions  concerning  the  propagation  of 
the  fundamental  mode. 

1.  What  is  the  size  of  the  focused  spot  obtainable  when 
a  sample  is  placed  in  the  beam? 

2.  What  happens  to  the  location  of  the  focused  spot 
when  a  sample  is  placed  in  the  beam? 

3.  How,  if  at  all,  are  the  spot  size  and  location 
affected  by  different  positions  of  the  air-sample 
interface? 

For  this  purpose  we  have  used  as  a  starting  point  an  equation 
derived  by  V.  Evtuhov*  from  the  theory  of  Boyd  and  Gordon.6 


W 


W2(r4+s)2  + 
- -  1 


a  2  2 

4s  r^ 


W. 

l 


(1) 


W2(ri+s)2  +  4s2r2 

5  Y~ 

W^(ri+s)  +  4sr^ 


(2) 
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where  r  is  the  radius  of  curvature  of  the  wave  front,  r^ 
is  the  initial  radius  of  curvature,  s  is  the  distance,  and  „ 
W  is  related  to  the  beam  radius  by  the  relation  W  =  (2tt/X)  w^ 
w  is  the  beam  radius  defined  as  the  radial  distance  for  2 
which  the  field  reaches  1/e  of  its  peak  value  (or  the  l/ez 
point  for  the  intensity) .  is  the  value  of  W  at  the 

origin  or  input.  If  we  have  a  gaussian  mode  propagating  in 
the  s  direction  starting  from  a  plane  through  s  =  0,  having 
an  initial  radius  of  curvature  rj_  and  beam  radius  w^,  eqs. 
(1)  and  (2)  give  the  value  of  W  and  r  at  any  distance 
s  from  the  input.  When  an  optical  element  such  as  a  lens, 
mirror,  or  dielectric  is  placed  in  the  beam,  we  use  the  values 
of  W  and  r  at  that  point  as  representing  new  initial  con¬ 
ditions  for  finding  further  values.  For  our  calculation, 
therefore,  we  start  at  the  laser  output  mirror  with  a  certain 
initial  beam  size  and  radius  of  curvature  (w^  =  0.1  cm,  r^  *  » 
and  allow  the  gaussian  wave  to  propagate  to  the  lens  a  dis¬ 
tance  s  (1  m  in  this  case)  away.  Using  (1)  and  (2)  we 
compute  the  value  of  r  and  W  at  the  input  to  the  lens.  At 
the  output  of  the  lens  W  is  unchanged,  whereas  r  takes  on 
a  new  value,  rou^  =  frj.n/  (f-r£n)  .  These  are  used  as  new  ini¬ 
tial  conditions  for  the  propagation  of  the  wave  after  the 
lens.  (Note  that  for  a  plane  wave  incident  on  the  lens 
(rin  =  “)  the  outgoing  radius  of  curvature  rput  =  -f.) 

When  a  material  of  dielectric  constant  n  is  placed  some¬ 
where  after  the  lens  it  is  again  necessary  to  compute  r  and 
W  at  this  point  and  use  these  as  new  initial  conditions.  In 
this  case  the  beam  radius  w  is  not  changed  at  the  boundary, 
but  the  radius  of  curvature  increases  by  a  factor  equal  to  the 
refractive  index  n  of  the  medium.  (More  precisely,  the 
radius  of  curvature  is  multiplied  by  the  ratio  of  refractive 
indices  of  the  two  media.) 

Let  us  define  the  following  unitless  quantities: 


Q  = 


A 


r 

r 


Equations  (1)  and  (2)  can  be  rewritten  in  the  following 
more  convenient  forms : 


n  =  (1  +  E)2  +  4  E 2  p  2 


(3) 


A  = 


(1  +  E) 2  +  4E2p2 
1  +  E  +  4  Ep2 


(4) 
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The  condition  for  minimum  beam  radius  can  be  found  by 
taking  .he  derivative  of  (3)  with  respect  to  l  or  setting 
the  denominator  of  (4)  equal  to  zero.  In  either  case,  one 
obtains 


^  (ft  min) 


-1 

2 

l+4p^ 


Further  substitution  into  (3)  gives 


ft  =  JlfiL. 

min  ..TT  * 

1+4  p 

Let  us  now  place  in  the  path  a  semi- inf inite  sample  of  re¬ 
fractive  index  n.  For  our  purposes  we  will  assume  that  the 
light  is  traveling  from  left  to  right.  The  sample  is  located 
such  that  the  air-sample  interface  is  found  at  an  arbitrary 
distance  b  from  the  lens.  To  the  right  of  boundary  b  (i.e., 
inside  the  material) ,  we  will  denote  the  parameters  as  primed, 
e.g.,  I',  p',  etc.  To  the  left,  the  variables  are  unprimed. 

At  the  boundary,  we  determine  new  initial  conditions 
for  the  further  propagation  of  the  wave.  The  radius  of  curva¬ 
ture  just  inside  the  boundary,  f£  =  nr^,  where  ru  is  the 
radius  of  curvature  just  outside  the  boundary.  Similarly, 
w£  =  wfa,  nX'  =  X;  hence,  W£  =  nW^. 

We  now  set  up  new  initial  conditions  at  the  boundary 
plane  b  under  two  conditions:  (1)  with  sample  present, 

(2)  without  sample. 

The  minimum  conditions  for  Case  1  are 


.  ) 

min 


S'  . 
min 
zrt 


-1 

- rj 

l+4pb 


ft  '  . 
min 


4Pv 


l+4p. 
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where 


p 


b 


p 


b  * 


The  conditions  for  Case  2  are  the  same,  except  that  we 
use  the  imprinted  values  corresponding  to  initial  conditions 
infinitesimally  to  the  left  of  the  boundary. 

We  see  that  the  value  of  ft  at  the  minimum  depends 
solely  on  the  value  of  p  which  is  the  same  on  either  side 
of  the  boundary.  Thus,  the  size  of  the  focused  spot  is  the 
same  regardless  of  whether  a  sample  is  present,  and  is  in¬ 
dependent  of  the  location  of  the  air-sample  boundary. 

Similarly,  the  values  for  Z  are  the  same  but  the 
distance  of  the  minimum  from  the  boundary  is  different  in 
the  two  cases. 

For  Case  1 


~nrb 


SmLn  1+4 pj 

Hence  »  nsm^n.  That  is,  the  distance  from  the  plane 

b  to  the  focus  is  n  times  farther  when  a  sample  is  present 
than  with  no  sample  present. 

We  summarize  the  results  of  the  preceding  considerations 
in  the  following  way.*  (1)  The  size  of  the  beam  waist  of  a 
focused  beam  is  unchanged  when  any  material  is  placed  in  the 


We  wish  to  acknowledge  the  comments  of  E.S.  Bliss,  who  pointed 
out  an  error  in  an  earlier  treatment  which  led  to  appreciably 
different  conclusions  from  those  we  now  reach. 


Sin 


”rb 


l+4p. 


For  Case  2 
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path.  (2)  The  location  of  the  beam  waist  is  moved  "down¬ 
stream”  n  times  farther  from  the  sample  boundary  than  it 
would  have  been  if  the  sample  had  not  been  present. 

H.  Plans  for  Next  Period 


Threshold  measurements  will  be  carried  out  under  a  variety 
of  different  conditions  as  soon  as  the  resonant  reflector  can 
be  installed.  The  parameters  which  will  be  varied  in  the  be¬ 
ginning  will  be  sample  length,  beam  size  (i.e.,  focal  length 
of  focusing  lens) ,  and  external  pumping.  The  last  may  be 
especially  important  if  photoconductivity  is  an  important 
factor  in  the  damage  mechanism.  Samples  from  different 
sources  will  be  subjected  to  test  as  soon  as  they  become  avail¬ 
able.  Among  these  will  be  TiO?  doped  ruby  and  sapphire  from 
Dj«va,  which  showed  a  considerably  higher  damage  threshold  than 
undoped  material  in  previous  work.7 
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II.  THEORETICAL  STUDIES  ON  OPTICAL  DAMAGE 


A.  Introduction  and  Summary  of  Results 


During  this  reporting  period  our  main  theoretical  efforts 
have  concentrated  on  the  effects  of  strong  optical  fields  on 
relatively  low  densities  lO^S  cm--*)  of  electrons  in  the 
conduction  band  of  a  transparent  polar  crystal  (with  special 
attention  to  ruby  and  sapphire  crystals)  .  Various  experi¬ 
menters  have  presented  evidence  that  electron  densities  of  at 
least  1()15  cm-2  are  photoexcited  in  ruby  and  sapphire  by 
optical  intensities  of  the  order  of  lO*0  W/cm2.  which  are 
just  below  damage  levels. *“,#  Previous  workers 1 • 1 1 “ 1 %  have 
suggested  that  at  a  certain  threshold  of  optical  pulse  power, 
these  electrons  become  hot  enough  to  ionize  their  surroundings 
and  produce  a  catastrophic  electron  avalanche  that  parallels 
the  electron  avalanche  process  long  known  and  studied  in  dc 
dielectric  breakdown.13  Our  own  studies  have  led  us  to  agree 
with  the  general  opinion  that  these  electrons  are  a  key  link 
in  the  optical-pulse  damage  process,  and  we  agree  with  all 
previous  workers  that  the  model  "polaron"  Hamiltonian  con¬ 
structed  by  Frdhlich  to  describe  the  electron- lattice  coupling 
should  give  at  least  a  good  qualitative  description  of  these 
nondegenerate  conduction  electrons.  However,  our  studies  of 
electron  heating  using  the  lowest  order  (in  coupling  and  opti¬ 
cal  intensity)  results  of  the  Frdhlich  model  do  not  indicate 
that  significant  electron  heating  above  ambient  lattice  tem¬ 
peratures  should  occur  at  the  optical  intensities  ('v  1010  W/cm2) 
in  question.  Furthermore,  we  find  that  if  these  lowest  order 
terms  give  an  adequate  description,  an  entirely  new  role  of 
the  electrons  in  damage  is  suggested:  the  electrons  cause  the 
optical  energy  to  be  passed  to  the  longitudinal-optical  (LO) 
phonon  modes  rapidly  enough  that  the  subsequent  decay  of  these 
LO  phonons  into  acoustic  phonons  (in  'v  10"^  sec)  sets  up 
acoustic  shock  waves  that  result  in  lattice  rupture.  We  have 
not  yet  studied  the  problem  of  whether  these  damage-mediating 
electrons  are  created  by  single-  or  multiphoton  excitation,  or 
possibly  by  photoemission  from  inclusions.  The  photoconduc¬ 
tivity  that  Hochuli*  observed  in  ruby  and  sapphire  was  linear 
in  the  optical  intensity  I,  so  that  it  occurred  presumably 
by  single-photon  excitation  from  impurity  levels  (a  presump¬ 
tion  supported  by  the  observed  temperature  dependence*).  How¬ 
ever,  the  photoconductivity  observed  in  ruby  by  Dneprovskii, 
et  al. ,*  at  much  higher  intensities  varied  as  I2.  In  either 
eventuality  it  would  appear  that  the  conduction  electrons  are 


being  drawn  from  impurity  levels.  Therefore,  if  the  new  pic¬ 
ture  of  the  damage  process  which  we  discuss  below  is  valid, 
there  is  real  hope  of  raising  the  damage  threshold  of  a  crys¬ 
tal  by  adding  impurities  which  bind  the  existing  impurity 
electrons  more  tightly  and  reduce  the  supply  of  photoelec¬ 
trons.  This  possibility  does  not  exist  if  electron  avalanche 
is  important,  because  avalanche  processes  occur  quite  inde¬ 
pendently  of  the  number  of  electrons  available  initially  to 
start  the  avalanche. 


B.  Estimate  of  the  Conditions  for  Crystal  Rupture 


If  the  role  of  the  electrons  is  as  we  have  suggested, 
their  presence  results  essentially  in  a  certain  thermal  (heat) 
energy  (which  we  calculate  in  the  next  section)  being  de¬ 
posited  in  the  lattice  by  an  optical  pulse.  To  estimate  an 
upper  limit  on  how  much  energy  is  required  to  rupture  the 
crystal,  we  use  the  observations  of  Belikova  and  Sviridenkov1 0 
of  the  damaging  shock  waves  that  formed  at  breakdown  around  a 
focal  region  ^  10-1  cm  long  and  ^  10“2  cm  wide.  By  measuring 
the  shock  velocity,  they  estimated  that  the  shock  energy 
('i  10_1  J)  was  essentially  the  same  as  the  energy  that  was 
observed  to  be  absorbed  and  scattered  from  the  0.3  -  0.5  J, 

10”8  sec  optical  pulse.  (They  estimated  that  this  energy  cor¬ 
responded  to  shock  front  pressures  of  the  order  of  the  yield 
stress  in  ruby,  105  psi.)  Therefore,  *3  may  estimate  that 
when  the  photoexcited  electrons  produce  an  optical  absorption 
coefficient  k  cm-1  of  order  unity,  damage  will  certainly 
result.  Other  observations  of  our  own  and  others1  have  indi¬ 
cated  that  in  fact  damage  may  occur  with  less  than  101  ab¬ 
sorption  from  the  beam;  we  believe  that  for  some  beam  geom¬ 
etries  (as  yet  uncertain) ,  as  little  absorption  as  k  ^  10“1 
to  10”2  is  sufficient  to  initiate  a  damaging  acoustic  shock. 

Unfortunately,  to  date  there  are  no  calculations  of  how 
the  damaging  shock  is  formed  from  the  temperature  and  pressure 
gradients  causod  by  the  deposited  energy,  so  that  we  have  no 
good  theoretical  corroboration  of  these  numbers.  However,  we 
have  been  able  to  estimate  from  standard  hydrodynamic  theory 
that  the  threshold  for  stimulated  thermal  Rayleigh  scattering 
from  a  10^°  W/cm2  beam  occurs  in  ruby  when  k  ^  1.  In  the 
next  section  we  show  from  polaron  theory  that  it  requires  on 
the  order  of  1016  electrons/cm3  in  a  ruby  or  sapphire  to 
produce  an  absorption  coefficient  of  order  10-3-  cm-1,  a  figure 
which  is  consistent  with  our  idea  that  the  photoelectrons  can 
damage  by  mediating  lattice  heating  without  ionizing  the  lattice. 
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c. 


Calculation  of  Heating  and  Absorption  by  Conduction 
Electrons 


Frohlich's  model  Hamiltonian  H  for  a  photoexcited  elec¬ 
tron  coupled  to  its  surrounding  lattice  and  in  an  oscillating 
x-directed  electric  field  E(t)  is 

H  ■  h  +  u  +  v  (5) 

where 

h  -  p2/2  +  £kVkak  (6) 


describes  the  uncoupled  electron  (of  momentum  p)  and  free 
longitudinal  polar  lattice  modes  of  frequency  w.  (whose 
annihilation  operators  are  a^) .  Here 

u  *  ^t,Ct,a.  exp(ik  •  x)  +  h.c.  (7) 

lv  T%  JC  ^ 

describes  the  electron-lattice  coupling;  and 

v  «  -xE  (t)  (8) 

describes  sufficiently  accurately  the  coupling  of  the  electron 
with  position  coordinate  x  to  the  optical  field.  We  employ 
polaron  units,  unless  otherwise  stated,  where  dS  ■  m*  ■  1, 

and  m*  is  the  electron's  (assumed  parabolic)  band  mass  for 
a  stationary  lattice.  Frohlich  derived  the  coupling  coeffi¬ 
cient  values  Ck  ■  23/4Trl/2al/2/k  where  the  dimensionless 

measure  a  of  the  electron  lattice  coupling  strength  is  de¬ 
fined  by 

a  =  (e_1  -  e"1)  (-fiai/Ryd)-1^2  (9) 

and  eQ,  t  are  the  dc  and  optical  dielectric  constants  of 
the  crystal.  The  longitudinal  optical  phonon  frequency  is  ap¬ 
proximated  by  the  constant  u  and  Ryd  is  the  Rydberg  con¬ 
stant  for  an  electron  of  mass  m*. 
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If  the  coupling  constant  a  is  not  large  competed  with 
unity,  calculated  quantities  may  be  developed  in  a  perturba¬ 
tion  series  in  a.  The  terms  of  order  a°  or  a  generally 
give  fairly  accurate  predictions,  even  for  a's  as  large  as 
unity.  For  ruby  or  sapphire,  t  ^  3,  eQ  ^  9,  but  the  effective 
band  mass  is  unknown.  If  we  take  the  band  mass  to  be  the  free 
electron  mass  and  take  iuo  ■  Boltzmann's  constant  k  times 
the  Debye  temperature  (980°K) ,  then  (5)  gives  a  ^  3,  which 
is  unfortunately  too  large  for  a  weak  coupling  approximation 
to  be  accurate.  To  examine  the  electron  energy  gain,  however, 
we  will  proceed  as  if  a  were  small  (as  did  all  previous 
theories  of  damage) ,  with  the  knowledge  that  methods  are 
available  to  calculate  this  more  accurately,1'  as  we  are  at¬ 
tempting  at  present. 

First,  to  determine  whether  the  electron  does  heat  up 
appreciably,  we  have  assumed  that  the  optical  field  was  turned 
on  at  t  *0  and  calculated  the  expected  kinetic  energy 
K  ■  1/2  <i2(t)>  of  the  electron  at  later  times  t  using  (1) . 
We  have  found  that 


2K  -  kT  +  V2(t)[l  +  01 (O)  +  02  (E2/E2)]  (10) 

where  T  is  the  local  lattice  temperature  and  V(t)  is  the 
expected  velocity  of  an  electron  initially  at  rest  and  driven 
by  the  electric  field  E(t)  in  the  presence  of  lattice  damp¬ 
ing.  The  corrections  Oi  and  O2  are  of  order  a  and 

(E/E  ) 2  where  E„  is  a  characteristic  electric  field 
o  o 

Eq  -  (w2/e)  (fun*/w)1/2  esu  (11) 

where  e  is  the  electric  charge  in  esu.  For  E  >  E0  the 
electron  develops  large  coherent  nonlinear  response.  This 
field  is  seen  to  be  on  the  order  of  10^  to  10°  V/cm  for  ruby, 
and  corresponds  to  optical  power  densities  of  the  order  of 
damage  threshold  densities.  (Our  calculations  of  these  non¬ 
linear  corrections  indicate  they  do  not  qualitatively  affect 
the  damage  process,  but  tend  to  decouple  the  electron  from 
the  lattice.  We  will  neglect  them  here.)  The  important  char¬ 
acteristic  of  all  the  terms  in  (6)  for  the  expected  kinetic 
energy  is  that,  after  a  brief  and  small  transient,  the  kinetic 
energy  becomes  essentially  constant  in  time  in  an  optical 
field  of  constant  amplitude.  This  is  in  contrast  to  the 


suggestions  of  all  previous  workers3'11-14  that  in  these  cir¬ 
cumstances  the  electron  might  continually  gain  energy  in  time 
until  it  had  enough  energy  ('v  10  eV)  to  ionize,  as  it  might 
indeed  gain  in  a  dc  electric  field.15  Because  no  previous 
worker  has  calculated  K(t)  and  examined  this  point,  we  can¬ 
not  point  to  any  calculational  error  that  might  have  led  them 
to  suppose  that  the  optical  damage  would  imitate  dc  dielec¬ 
tric  breakdown  in  that  very  hot  electrons  would  be  created. 

It  is  easy  to  estimate  that,  since  (as  we  will  show)  the 
electron-to-lattice  damping  rate  is  small  compared  with  the 
optical  frequency  v,  V(t)  E(t)/v2  so  that  <  V2>  is  of  the 
order  of  10"5  eV  when  E  is  near  the  damage  level.  There¬ 
fore,  to  the  lowest  orders  of  approximation  the  average  elec¬ 
tron  kinetic  energy  is  always  close  to  the  thermal  energy 
('i  1/40  to  1/4  eV)  that  it  would  have  in  equilibrium  with  the 
lattice  in  its  vicinity.  Finding  this  result,  we  proceeded  to 
calculate,  using  (5) ,  the  expected  rate  of  increase  in  the 
total  phonon  energy  to  lowest  order  in  a.  We  found,  as  ex¬ 
pected,  that  the  phonon  energy  does  increase  proportionally 
with  time  (in  a  constant  optical  field)  with  a  rate  that 
matches  the  rate  calculated  by  others14  at  which  the  electric 
field  does  work  on  the  electron. 

This  confirms  the  idea  that  the  electrons  nearly  instant¬ 
ly  pass  all  the  energy  that  they  receive  from  the  field  on  to 
the  optical  phonons.  We  have  not  yet  been  able  to  confirm 
that  the  field-to-electron  energy  transfer  rate  matches  the 
electron-to-lattice  rate  at  arbitrary  coupling,  temperature, 
and  optical  intensity,  but  the  structure  of  the  general  equa¬ 
tion  for  K  indicates  that  it  is  indeed  constant  in  a  con¬ 
stant  amplitude  field,  and  thus  we  now  believe  that  the  crea¬ 
tion  of  hot  electrons  when  they  are  coupled  only  to  a  cold 
lattice  is  improbable.  However,  a  possible  complication  that 
we  will  study  further  is  that  a  very  few  electrons  might  be 
created  in  the  high  energy  wings  of  the  electron  distribution. 

To  obtain  an  estimate  of  the  expected  rate  y  at  which 
the  optical  field  of  amplitude  E  oscillating  at  a  frequency 
v  does  work  on  an  electron,  we  may  use  the  polaron  response 
theory  derived  by  Feynman,  et_  al . , 1 6  using  Frohlich's 
Hamiltonian 


Y  -  \  E2  lm[v2  -  x (v) )-1  (12) 
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and 


cos(vu)  co8u  du _ 

[u2  +  a2  -  b  cos  (vu)  ]  3//2 
o 

(13) 

the  phonon  frequency  w 
is  1  (making  v  ^  17  ?.t  the  ruby  laser  frequency)  ,  time  is 
expressed  in  units  w“l,  energy  in  units  of  -hw,  etc.  The 
parameter  $  =  hoi/kT  is  ^  3.3  for  ruby  and  sapphire  at  room 

temperature.  The  parameters  v  and  w  are  chosen  from  a 
variational  principle  which  gives  v  =  3.4,  w  =  2.55  for  a  *  3 
and  in  the  limit  of  low  temperatures16;  a2  =  B2/4  +  R8  coth 
(8v/2) ;  R  =  (v2  -  w2 ) / (w2v) ;  and  b  =  R$/sinh(Bv/2)  .  With  the 
above  parameters  R  =  0.229,  a  =  1.866,  and  b  =  5.53  x  10”3. 

The  integrand  in  (13)  can  be  expanded  in  powers  of  b 
and  each  term  integrated  to  obtain  the  value  of  y.  Unfortu¬ 
nately,  the  series  converges  rather  slowly  for  the  high  fre¬ 
quencies  of  interest,  reflecting  the  fact  that  the  oscillating 
electron  can  lose  energy  by  processes  in  which  1,  2,  or  many 
phonons  are  omitted.  We  are  in  the  process  of  evaluating  the 
integral  numerically  for  various  parameters.  A  first  estimate 
of  (13)  using  the  above  parameters  shows  Imy  ^  40  ,  a  value 
sufficiently  accurate  for  our  present  qualitative  estimates. 

The  optical  beam  absorption  coefficient  k  is  simply  re¬ 
lated  to  x  by 


Imx(v)  = 


2a7T~1/283/2sinh(gv/2)  /v\3 


3  sinh(8/2) 


GO 


We  have  redefined  our  units  so  that 


K/ko  =  (u)p/e)  Im  [y2  ~  X)"1  (14) 

where  (*>„  is  the  electron  plasma  frequency  in  our  units  and 
equals  (47rNe2/m*) l/2  in  cgs  units,  where  N  is  the  number 
density  of  electrons.  To  obtain  k  ^  0.1  cm"l  when  kQ,  the 
optical  wavevector  length  in  the  crystal  is  103  cm”l,  we 
would  require,  according  to  our  above  estimates,  that  u>p  be 
on  the  order  of  1/30  of  the  phonon  frequency.  This  corresponds 
to  an  electron  density  of  ^  10^^  cm-3  if  the  effective  mass  is 
unity.  It  is  important  to  note  that  with  ImX/  there  is  a  real 
part  of  x  which  will  contribute  to  a  nonlinear  index  and 
which  also  has  the  effect  of  altering  the  effective  mass  of 
the  electron,  although  probably  not  by  more  than  a  factor  of  2. 
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D. 


Plans  for  Next  Period 


To  examine  the  validity  of  our  new  picture  of  the  role  of 
photoexcited  electrons  in  crystal  damage,  we  plan  to  (a)  to 
try  to  establish  a  firm  dynamical  acoustic  theory  for  the  for¬ 
mation  of  shock  waves,  given  the  thermal  energy  profile  de¬ 
posited  by  the  optical  beam;  (b)  to  evaluate  precisely  the 
predictions  of  Frohlich's  polaron  model  for  the  damping  of 
the  driven  optical  motions  of  the  electrons;  (c)  to  evaluate 
the  tendency  of  the  electron  to  gain  energy  in  the  form  of 
incoherent  motions  under  the  influence  of  the  optical  beam; 

(d)  to  examine  further  the  various  mechanisms  of  photoexcita¬ 
tion;  (e)  to  calculate  the  nonlinear  index  of  refraction 
seen  by  the  optical  beam  because  of  the  presence  of  conduction 
electrons;  (f)  to  calculate  the  nonlinear  effects  on  the  ab¬ 
sorption  rate  that  occur  at  high  fields;  and  (g)  to  continue 
to  search  for  possible  processes  in  which  these  electrons  are 
not  important  to  damage  at  all. 
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with  the  general  opinion  that  these  electrons  are  a  key  link 
in  the  optical  pulse  damage  process.  However,  the  present  work 
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temperatures  does  not  occur  at  the  optical  intensities  in  ques¬ 
tion.  An  entirely  new  role  of  the  electrons  in  laser  induced 
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